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Abstract. In a recent Letter, Cowsik, Ratnam and Bhat- 
tacharjee (1996a) have built a dynamically self-consistent 
spatial distribution of particles of galactic dark matter. 
They have come up with the rather unorthodox conclu- 
sion that the mean velocity dispersion of dark matter 
particles "should be 600km s _1 or larger". Their letter 
triggered immediate comments (Evans 1996, Gates et al. 
1996). Here we find, as did Cowsik et al. (1996a), that 
models of the dark matter halo can be made consistent 
with velocity dispersion much larger than that expected 
from a simple application of the virial theorem in the so- 
lar neighbourhood. But in contrast to their conclusions, we 
show that using their model, we also obtain solutions with 
smaller velocity dispersion down to ~ 270 km s _1 . These 
more orthodox dispersions arise because of less constrain- 
ing boundary conditions for the central density but do not 
rely on indirect or model dependent measurements of the 
large scale distance behaviour of the rotation curve. 
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We developed a numerical algorithm following the pre- 
scription presented by Cowsik et al. (1996a) and we obtain 
a rotation curve with the same general features. This re- 
sult provides an independent validation of Cowsik et al.'s 
numerical work. The distribution of the baryonic matter 
in our Galactic model corresponds to a double-exponential 
disk (with a scale length of 3.5 kpc, a scale height of 300 pc 
and a surface density, So, of 80 M Q pc~ 3 at the solar ra- 
dius i? = 8.5 kpc). The bulge 's model is a Hubble pro- 
file with a core radius a = 103 pc and a central density 
p c = 343 A/ Q pc- 3 . 

We check that the velocity curves of these models do 
tend asymptotically to (2/3) 1 / 2 (w 2 ) r / I ^ = ©oo- We de- 
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termine velocity curves for three models with different 
(v 2 ) 1 ^ = 350, 450 and 600 km s" 1 (Fig. 1) (the fixed size 
of our rectangular grid does not allow us to resolve simul- 
taneously with sufficient accuracy both the velocity curve 
to galactic radius of 40 kpc and close to the galactic center 
below 2 kpc which we have therefore ommitted) . 

For each dispersion we estimate values of pdm{Ro = 
8.5 kpc, 0) in order to obtain a flat rotation curve in the 
range R = 10 — 30 kpc. We note a dip at about 20 — 25 kpc, 
which increases with the velocity dispersion (too high a 
velocity dispersion prevents a nearly flat rotation curve). 
The velocity curve maximum near R = 7 kpc is produced 
by the exponential disk of visible matter. At much larger 
galactic radius R > 40 kpc the dark matter contribution 
dominates the velocity curve which tends towards ©oo. 
For instance in Fig. 1, the lowest curve (3) that corre- 
sponds to the largest {v 2 }^^ reaches 490 km s _1 at large 
R while curves (2) and (1) only reach 367 and 286 km s _1 . 
Increasing pdm(Ro,0) increases the mass contribution of 
the dark halo in the central region and convergence to- 
wards ©oo is reached more rapidly. Yet, in this instance, 
the velocity curve does not remain flat below 40 kpc. 
The most unusual feature of the models proposed by 
Cowsik et al. is that the central dark matter density is so 
low that the visible matter dominates the velocity curve at 
small galactic radius and allows for a galactic population 
with an extension stopping at 30 — 40 kpc while its veloc- 
ity dispersion is (2/3) 1 / 2 220 = 180 km s _1 (very similar to 
that of the Globular Cluster population for instance). The 
dark matter population having a larger (v 2 )^ dominates 
fully after 30 kpc and as shown in Cowsik's reply (1996bc) 
this dark matter component may stop at 100 kpc without 
significant effect at small galactic radius. 

Cowsik's results are in conflict with available data 
only in the range 20 — 40 kpc. In this range, constraints 
arise from the dynamics of distant objects such as globu- 
lar clusters. These give access to the galactic mass inside 
R ~ 40 kpc, and suggest that its distribution is consistent 
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Fig. 1. Flat rotation curve of the Galaxy (f 2 )^ ~ 350 (1), 450 (2) and 600 km s _1 (3). The observational data and their error 
bars are those given by Cowisk et al : from Fich et al (1989), below 2kpc from Burton & Gordon (1978) and for R > 17kpc 
from Fich & Tremaine (1991). The corresponding \ 2 take respectively the values of 14, 18 and 22km s _1 for models (1), (2) 
and (3) 



with a flat rotation curve up to that distance (Dauphole 
& Colin 1996). In this reference, input data are 6D (r,v) 
positions and velocities for a critical subset of globular 
clusters, giving for the first time strong constraints on the 
rotation curve in the analysed range of distances. Their 
most probable fit corresponds to a flat rotation curve, 
but no upper limits to a rising rotation curve are given. 
Of course, such results remain model dependent and do 
not exclude a rising (or decreasing) rotation curve beyond 
40kpc. 

We remark that the data used by Cowsik et al. (1996a) 
correspond to a rotation curve determined with the as- 
sumption that the solar galactic radius is Rq = 8.5 kpc (a 
value recommended by the IAU - see Kerr & Lynden-Bcll 
1986), used to "facilitate intercomparison of the work of 
the different authors" (Fich & Tremaine 1991). Smaller 
values (Reid 1993) are now generally accepted and lead 
to a flat or locally decreasing rotation curve (Fich & 
Tremaine 1991). 

Our models do not take as data input the mass de- 
termination from 17 satellites of the Galaxy more dis- 
tant than 50 kpc (Kochaneck, 1996). Only 3 have complete 
(r, v) data and mass estimate depends on complementary 
hypotheses like velocity isotropy or radial gradient of the 
density. Similarly we do not incorporate the timing argu- 
ment as means to estimate the total mass of the Galaxy 
since as argued by Kochaneck (1996) " The classical Local 
Group timing model of Kahn and Woltjer 1959 assumes 



that the orbits are radial and provides lower bounds on 
the mass." 

Within the framework of our isothermal model, we find 
that an upper limit of 600 km s _1 is consistent with a flat 
rotation curve in the range R = 10 — 40 kpc. A range of 
possible galactic dark matter velocity dispersions from 270 
to 600 km s _1 are shown to be also consistent. We remark 
that the velocity dispersion can be much larger beyond 
40 kpc without significant measurable effects on the rota- 
tion curve below this radius. If the velocity curve of the 
Galaxy is only well defined over such a small distance in- 
terval, it is difficult to model it with a unique isothermal 
model. One should therefore be careful in estimating ve- 
locity dispersions from a crude application of the virial 
theorem, as mentioned by Cowsik et al. (1996a). 

We find a large range of possible (Odm with dis- 
persion as low as 270 km s _1 if we consider the ro- 
tation curve used by Cowsik et al., and from 270 to 
600 km s -1 if we consider a flat rotation curve up to 
40 kpc. The discrepancies between Cowsik's et al. con- 
clusion and ours arises from the fact that we did not 
fix the dark matter density at the origin pd.m. (0,0). 
Cowsik et al. set pd.m. (0,0) = lGev/cm 3 constant, ar- 
guing that the dynamical measure of the mass density 
in the solar neighbourhood is ~ 0.3 Gev/cm 3 . This value 
of Pd.m.(Ro,0) ~ 0.3 Gev/cm 3 = 0.008 M e /pc 3 is very 
different from the values (0.05 — 0.1 M Q /pc 3 ) obtained in 



O. Bienayme & C. Pichon: Comment on the dispersion-velocity of galactic dark matter 



3 



references (Oort 1960, Bahcall 1984) cited by Cowsik et al. 
(1996a). The exact estimate of the local dynamical mass 
density pDyn.{Ro, 0) remains a very controversial subject, 
though it is now determined with smaller errors (see re- 
view by Creze 1991 & Kuijken 1995). Models plotted in 
Fig. 1 have small pd.m. (Ro, 0) compatible with this last 
determinations of poyn. (Ro, 0); these models are in rough 
agreement with the data used by Cowsik in the range 
R < 20kpc and have much smaller (w 2 )^ values than 
those obtained by Cowsik et al. (1996a). 

Cowsik et al. have pointed out the important fact that 
models of dark matter halo can be made consistent with 
velocity dispersion much larger than that expected from a 
simple application of the virial theorem in the solar neigh- 
bourhood. Their model should be applied to existing 6D 
data from globular clusters in order to obtain a realis- 
tic upper limit on (w 2 )J/m- Here wc have shown that, in 
contrast to Cowsik et al. assertion, models with small or- 
thodox (w 2 )J/m ~ 270 km s _1 are also in agreement with 
their own observational constraints. Our conclusions sug- 
gest that the efforts on the part of experimentalists should 
not be directed towards the search for very hypothetical 
high velocity dispersion dark matter particles. These effort 
should on the contrary be directed towards better astro- 
metric observations in order to reach valuable absolute 
proper motions for the most distant galactic satellites and 
globular clusters that may be achieved for instance with 
high resolution observations with HST or adaptativc op- 
tics. 
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